Effect of vibrations on the pre-edge features of x-ray absorption spectra 
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The influence of atomic vibrations on x-ray absorption near edge structure (XANES) is calculated 
by assuming that vibrational energies are small with respect to the instrumental resolution. The 
resulting expression shows that, at the Tf-edge, vibrations enable electric dipole transitions to 3s 
and 3d final states. The theory is applied to the K-edge of Al in Q-AI2O3 and of Ti in Ti02 rutile 
and compared with experiment. At the Al if-edge, sizeable transitions towards 3s final states are 
obtained, leading to a clear improvement of the agreement with experimental spectra. At the Ti 
_K"-edge, electric dipole transitions towards 3d final states explain the temperature dependence of 
the pre-edge features. 

PACS numbers: 78.70.Dm, 65.40.-b 



Vibronic coupling describes the interaction between 
electrons and atomic motions. It plays a prominent role 
in optical spectroscopy where it is the source of the color 
of many pigments and gemstones For instance, the 
red color of Black Prince's ruby is due to "d-d" tran- 
sitions of chromium impurities in a spinel crystal. But 
these transitions are forbidden because chromium occu- 
pies an inversion center of the spinel lattice. They be- 
come allowed when vibrations break inversion symmetry. 

In the x-ray range, vibrations far from the edge 
are taken into account through a Debye- Waller factor 
e -2/c 2 cr 2 pj if ^he validity of this factor is assumed to 
extend to the near-edge region, where k ~ 0, then vibra- 
tions seem negligible in XANES spectra. 

However, three arguments indicate that vibronic cou- 
pling can be sizeable in the XANES region: (i) Vibronic 
coupling was detected by x-ray resonant scattering exper- 
iments at the Ge K-edge 0,H|; (ii) Some XANES peaks 
seem to be due to forbidden transitions to 3s states, a 
prominent example being the Al K-edge in minerals [5|; 
(iii) A temperature dependence of the pre-edge structure 
was observed at the Ti K-edge in SrTi0 3 @ and Ti0 2 Q. 

In the optical range, the effect of vibrations is usually 
taken into account through the Franck-Condon factors. 
In the x-ray range, Fujikawa and coll. showed in a series 
of papers of increasing sophistication 0, 0, that the 
effect of the Franck-Condon factors can be represented 
by the convolution of the "phonon-less" x-ray absorp- 
tion spectrum with the phonon spectral function. Such 
a convolution leads to a broadening of the peaks with in- 
creasing temperature but this effect is hardly observable 
in the pre-edge region. 

Moreover, in the x-ray range it was shown that the 
large core-electron-phonon coupling of the Is core hole 
of carbon in diamond induces a strong lattice distortion 
and significant anharmonic contributions 11 1 . 

The key observation is that all these effects can be eas- 
ily taken into account if the vibrational energies are small 
with respect to the XANES spectral resolution (core hole 



lifetime -I- instrumental resolution) . This condition is cer- 
tainly not satisfied at the C K-edge [ll| but it becomes 
reasonable at the Al and Ti K-edges. In that range the 
XANES resolution is around one eV whereas the energy 
of vibrational modes is of the order of a few hundredths 
of eV although, of course, several phonons can be simul- 
taneously present. 

In this paper, we first use this observation to derive 
a manageable expression for the vibronically-coupled x- 
ray absorption spectra. Then, we apply the so-called 
crude Born-Oppenheimer approximation to further sim- 
plify this expression, so that only the core-hole motion 
remains. The resulting equation is compared with ex- 
periment in two different cases. At the Al K-edge in a- 
AI2O3 (corundum), vibrations induce transitions to 3s fi- 
nal states. These (Is — > 3s) monopole transitions explain 
a pre-edge peak that is completely absent from standard 
calculations. At the Ti K-edge in Ti02, Is — > 3d transi- 
tions are induced by vibrations. This explains why only 
the first two pre-edge peaks grow with temperature. 

XANES formula within the Born-Oppenheimer frame- 
work. According to the Born-Oppenheimer approxima- 
tion [HI, EH , the wavefunction of the system of electrons 
and nuclei can be written as the product x^(Q)^? n (r; Q)i 
where r is the electronic variable and Q = (Qi,...,Qjv) 
collectively denotes the position vectors of the N nuclei of 
the system. The electronic wavefunction ^> n (r;Q), with 
energy e„(Q), describes the state of the electrons in a 
potential where the nuclei are fixed at position Q. The 
ground state corresponds to n = 0. The origin of the nu- 
clear variables is chosen so that Q = is the equilibrium 
position, i.e. eo(0) is the minimum of eo(Q). 

For each n, the vibrational wavefunctions Xn(Q) are 
the orthonormal solutions of the Schrodinger equation 

(H kin (Q) + e„(Q)) X ^(Q) = EixtiQ). 

The total energy of the electrons + nuclei system is E^. 
Within the Born-Oppenheimer approximation, the x-ray 
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absorption cross-section is 

a{uj) = 4^ 2 a ^^| /'dQdrx}(Q)> / (r,Q)*£T 

Si 

xxo(Q)^o(r,Q)| 5{Ej-E -hw), 

where ao is the fine structure constant, fku the energy of 
the incident x-rays and e their polarization vector. The 
core-hole lifetime and the instrumental resolution can be 
represented by the convolution of the absorption cross- 
section with a Lorentzian function (T/7r) / (ui 2 +T 2 ). This 
gives us 

a» = 47ra ^| / dQdr X }(Q)V/(r, Q)*£ • r 



Si 



; Xo (Q)7A (r,Q) 



(Ej -E - fiw) 2 + 7 2 



where 7 = fiT. 

The energy Ej, can be written as the sum of the elec- 
tronic energy at equilibrium position e f and a vibrational 
energy Ej = e/ + E^ h . When 7 is much larger than 
the vibrational energy we can neglect the contribution 
of E^( h and sum over the vibrational states X/(Q)- The 
completeness relation gives us 

E4(Q)V/(q') = 5(q-q')- 



Therefore, 

cr 7 (w) = 4ira a / dQ^I / dr^/(r, Q)*e • r 
J / 

(e/ - e )7 



;Xo(Q)^o(r,Q) 



(e/ - e - fiw) 2 + 7 2 



•(I) 



Note that we derived this result without making the har- 
monic approximation. Therefore, the possible anhar- 
monic behavior due to the core hole [111] is taken into 
account. Note also that the final and initial energies e/ 
and eo do not depend on Q. Therefore, eq. (TTJ) is not 
the average of standard XANES spectra over various nu- 
clear positions Q. In other words, we have here a way to 
distinguish thermal disorder from static disorder due to 
impurities and structural defects. 

Now we make a different approximation for the initial 
and final electronic states. For a K-edge, the Is core level 
wavefunction is highly localized around the nucleus and 
it weakly depends on the surrounding atoms. Therefore, 
we can approximate ipo{r, Q) by <fio(r— Q a ), where Q a is 
the position vector of the absorbing atom and where </>o 
is the Is wavefunction of the absorbing atom at equilib- 
rium position. For the final electronic states ^/(r, Q) in 
the presence of a core hole, the variation of the nuclear 



coordinates Q in eq. |T|) is ruled by the vibrational wave- 
function xo of the initial state, which is expected to be 
rather smooth. Therefore, we make the standard crude 
Born-Oppenheimer approximation, according to which 
the electronic wavefunction does not significantly vary 
with Q for small vibrational motions. In other words, 
V>/(r, Q) ~ 0/(r), where 4>f(r) — ipf(r, 0). This gives us 



a j (to) 



Airaohu; J dQ| X o(Q)| 2 ^| J dr^/(r) 



e ■ r 



x0 o (r- Q Q ) 



/ 



(ef - e - hu) 2 + 7 2 



When the crude Born-Oppenheimer approximation is not 
valid, it is possible to Taylor-expand ipf{r, Q) as a func- 



tion of Q [14j. The integral over electronic variables de- 



pends only on the position of the absorbing atom, from 
now on denoted by R. Therefore, we can integrate over 
the other nuclear variables and the expression becomes 



4na hiu J dRp(R)^| J dr%(r)' 



e ■ r 



X(f> (r - R) 



/ 



(e/ - eo - ftw) 2 + 7 2 ' 



(2) 



where p(R) = / dQ| X o(Q)| 2 <5(Qa - R). Within the har- 
monic approximation, the core displacement distribution 
has the form 11511. 



p{R) = exp ( - 



R U 



R 



where U is the thermal parameter matrix [16| that is 
measured in x-ray or neutron scattering experiments. 

Calculation of the matrix element. For a hydrogenoid 
atom, the Is core-hole radial wavefunction is propor- 
tional to e~ ar , where a = Z/a , Z is the atomic number 
and ao the Bohr radius. The Is wavefunction (f>o(r) of a 
true atom is close to that of a hydrogenoid one and can 
be written as a fast converging linear combination of ex- 
ponentials. Thus, </>o(r — R) becomes a sum of shifted ex- 
ponentials that can be described by the Barnett-Coulson 
expansion 17 1 



e -a|r-R| = J2(2n + l)P n (r-R)c n (r,R), 

n 

with P n a Legendre polynomial and 

c n {r,R) = ^=(r < / r ' l+1/2 (ar<) J ft: n+1/ 2(ar > ) 

\/rR K 

+r> 4+1/2 (ar< )K' n+l/2 (ar> ) J , 



where r< (r>, resp.) is the smaller (larger, resp.) of r 
and R, I v {z) and K v (z) are the modified Bessel functions 
and I' v [z) and K' v (z) their derivatives with respect to 
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z. For notational convenience, we consider that the core 
wavefunction can be represented by a single exponential 
<j> (r) = Ce~ ar . 

To calculate the matrix element, we expand the fi- 
nal state wavefunction over spherical harmonics tpf(v) = 
J2e m fem(r)Y™(r). The matrix element over the elec- 
tronic variable is 

f di%(r)*e.r<Mr,R) = E X "( R )' 



Cm 



with 



oo . 

(2n + l)c n (r, R)P n (r ■ R). 
Standard angular momentum recoupling leads to 
(4tt) 2 C 



£ Jr 3 drft m (r)c n (r,R) 



n=\£±l\ 



J2(-l) m Yr X (e)Y*-™(R)C 



Ira 

lXnm — X 



, (3) 
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FIG. 1: Experimental |23j ] and calculated Al if-edge isotropic 
spectrum of corundum at 300 K. (7 C xp and Urn are the exper- 
imental and fitted thermal parameters (see text). 



where C*" are Gaunt coefficients. Equation ([3]) shows 
that all values of the final state angular momentum £ are 
now allowed. The core-hole wavefunction is still spher- 
ical, but with respect to a shifted centrum. Thus, with 
respect to the original spectrum, it is a sum over all an- 
gular momenta given by the Barnett-Coulson expansion. 
Thus, all final states angular momenta are available in 
spite of the fact that only electric dipole transitions are 
allowed. In particular, vibrations allow for dipole transi- 
tions to the 3s and 3d final states at the -FT-edge. 

General features of vibrational transitions. The fore- 
going approach enables us to draw some general conclu- 
sions concerning the effect of vibrations on XAS pre-edge 
structure. This effect is measurable if the density of non-p 
states of the system in the final state (i.e. in the presence 
of a core hole) is large and well localized near the Fermi 
energy (vibrational transitions towards p-states would be 
masked by the allowed vibrationless transitions) . For ex- 
ample, just above the Fermi level, many aluminum or 
silicon compounds have a strong density of 3s states and 
many transition metal compounds have a large density of 
3d states. In the first case, vibrational transitions appear 
as monopole Is — > 3s transitions, which are completely 
excluded with electromagnetic transitions. In the second 
case, vibrational transitions superimpose upon electric 
quadrupole Is — ► 3d transitions. Thus, vibrations induce 
transitions at specific energies in the pre-edge but hardly 
modify the rest of the XANES spectrum. 

A temperature dependence of vibrational transitions 
is expected if U varies with temperature. This occurs 
between K and room temperature if the sample has soft 
modes (i.e. low energy phonons). Otherwise, vibrational 



transitions are only due to the zero-point motion of the 
nuclei. In other words, vibrational transitions are then a 
consequence of the fact that, even at K, nuclei are not 
localized at a single point. 

We test these conclusions with two examples: the Al 
-KT-edge in corundum (where vibrational transitions to 3s 
states are expected) and the Ti K-edge in rutile (where 
vibrational transitions to 3d states and a temperature 
dependence are expected). 

The Al K-edge in corundum. The X-ray absorption 
cross section within the crude Born-Oppcnheimer ap- 
proximation has been implemented in the XSpectra pack- 
age [HI of the Quantum-espresso suite of codes [l^l • To 
calculate the integral over R in eq. @, it is found suf- 
ficient to compute the integral over a cube of size 2A, 
where A is the largest eigenvalue of the matrix U. This 
cube is cut into 27 smaller cubes where the integral is 
carried out using eq. (25.4.68) of Ref. |2fJ. The technical 
details of the self-consistent calculation are the same as 
in Refs. [13, H|. 

Figure Q] shows the result for experimental mean dis- 
placements (Ti = cr 2 =0.048 A, cr 3 =0.049 A 0j. The 
mean displacements at are the square root of the eigen- 
values of U and they have a more direct physical meaning 
than U. The vibration transitions are observed exactly 
at the position of the pre-edge peak which is absent from 
the calculation without vibrations. However, the vibra- 
tional transitions are overestimated and a better agree- 
ment is obtained when setting the mean displacements 
to 0.026 A. Therefore, computed vibrational transitions 
show up at the right position but with a too large inten- 
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FIG. 2: Experimental [3] and calculated Ti K pre-edge spec- 
trum of rutile. El and E2 are electric dipole and quadrupole 
transitions. 



and room temperatures, respectively. These values are 
more than ten times smaller than the experimental val- 
ues [Hj]. Moreover, when the thermal factor is adjusted 
to be consistent with the temperature variation of the 
second peak, then the temperature variation of the first 
peak is overestimated. This can probably be attributed 
to the crude Born-Oppenheimer approximation. Despite 
these drawbacks, several aspects of the experimental vi- 
brational transitions are correctly reproduced: (i) only 
the first and the second peaks exhibit a temperature de- 
pendence, the rest of the spectrum is not modified; (ii) 
the peaks do not shift and do not broaden; (iii) the peaks 
increase with temperature. 

This work was performed using HPC resources from 
GENCI grant 2009-2015 and 1202. 



sity. A part of this discrepancy might be due to the fact 
that the experimental thermal parameters include some 
amount of static disorder, but most of it probably comes 
from the crude Born-Oppcnhcimcr approximation. The 
main vibrational transitions come from the relative dis- 
placement of the aluminum atom with respect to the six 
oxygen first neighbors. The vibrational modes involving 
an overall motion of the aluminum atom with its oxy- 
gen octahedron contribute to the thermal parameter U 
although they do not contribute to the vibrational transi- 
tions. Because of the crude Born-Oppenheimer approx- 
imation, the vibrations are summarized in the thermal 
parameter, which overestimates the real effect of vibra- 
tions. A similar phenomenon occurs with the EXAFS 
Debye- Waller factor which is different from the one de- 
termined by x-ray diffraction because only relative dis- 
placements must be taken into account. This point will 
now be confirmed with the Ti K-edge absorption in ru- 
tile. 

The Ti K-edge of rutile. With the Ti K-edge of rutile, 
we test the limit of the crude Born-Oppenheimer approx- 
imation. Indeed, this approximation assumes that the fi- 
nal state wavefunction does not change when the crystal 
vibrates. This might be reasonable for the Al 3s states 
because they are poorly localized and overlap the oxygen 
2p orbitals, but this is not true for the 3d states of Ti in 
rutile, which are localized near the Ti nucleus. 

Figure [2] shows the experimental and theoretical Ti 
K-edge spectrum of rutile with the polarization parallel 
to the c axis and at two temperatures. A similar re- 
sult is obtained when e is perpendicular to the c axis. 
The mean displacements o~\ , 02 , 03 used in the calcula- 
tions are, in A, 0.0028, 0.0027, 0.0023 for Ul (8 K) and 
0.0043, 0.0042, 0.0035 for U2 (250 K). They have been 
chosen to approximate the intensity of peak A2 at low 
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